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We study the transient behavior in coupled dissipative dynamical systems based on the linear 
analysis around the steady state. We find that the transient time is minimized at a specific set of 
system parameters and show that at this parameter set, two eigenvalues and two eigenvectors of 
Jacobian matrix coalesce at the same time, this degenerate point is called the exceptional point. For 
the case of coupled limit cycle oscillators, we investigate the transient behavior into the amplitude 
death state, and clarify that the exceptional point is associated with a critical point of frequency 
locking, as well as the transition of the envelope oscillation. 

PACS numbers: 05.45.Xt, 02.10.Ud 


I. INTRODUCTION 

In the eigenvalue problem of a non-Hermitian matrix, 
an exceptional point (EP) is a square-root branch point 
on a two-dimensional parameter space, at which not only 
^envalues but also the associated eigenvectors coalesce 
[il, 0 . The peculiar feature related to the EP is the ex¬ 
change of eigenvalues and eigenvectors after a parame¬ 
ter variation encircling the EP once, of which topological 
structure is same as that of Mobius strip Q . The EPs and 
relating interesting phenomena have mainly been studied 
in open quantum systems described by non-Hermitian 
Hamiltonians such as atomic spectra in fields Hi , mi¬ 
crowave cavity experiments p, Q , chaotic optical micro¬ 
cavities [H, PT-symmetric quantum systems IM , and 
so on. Besides the open quantum systems, the EPs are 
also observed in coupled driven damped oscillators real¬ 
ized by electric circuits, which are purely classical sys¬ 
tems Bin. 

The amplitude death (AD) is the complete suppression 
of oscillations of the entire system when the nonlinear dy¬ 
namical systems are coupled H- The AD has been ob¬ 
served in many coupled dynamical systems and the AD is 
achieved by various types of coupling interaction, i.e., the 
diffusive coupling in mismatched oscillators llbl - llq . de¬ 
layed couplin g [l9l - [^ . conjugate couplin g 1 ^ . dynam¬ 
ical coupling |25l|. nonlinear coupling l27j. etc. The 
AD has also been studied in networks of coupled oscilla¬ 
tors [13, HH and variety topologies such as a ring [29|, , 

small world [3l|, and scale free networks [32|. Recently, 
the suppressions of oscillations are strictly classified into 
amplitude death and oscillation death, where the asymp¬ 
totic steady state is homogeneous and inhomogeneous, 
respectively. [ 13 , [33| 

In this paper, we study the transient behaviors of cou¬ 
pled dissipative dynamical systems based on the linear 
analysis around the steady state. We find that the sys- 
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terns show the largest damping rate at an EP, which 
comes from the intrinsic feature of a square-root branch 
point. Eor the case of coupled limit cycle oscillators, the 
transient behavior into the amplitude death state is stud¬ 
ied. We demonstrate that the EP is associated with a 
critical point of frequency locking, as well as the transi¬ 
tion of the envelope oscillation. 

This paper is organized as follows. In Sec. H, we show 
the occurrence of EP in coupled damped oscillators and 
discuss the damping behavior around the EP in a peda¬ 
gogical way. In Sec. HI, we present the transient behav¬ 
ior into the AD in coupled limit cycle oscillators, and it 
is explained based on the existence of an EP. Einally, we 
summarize our results in Sec. IV. 


II. EXCEPTIONAL POINT IN COUPLED 
DAMPED OSCILLATORS 

We consider the coupled damped oscillators, 

^ 1 + 71 ^ 1 +^ 1^1 = -kx2, 

^2 + 72^2 +^ 2^2 = -kXi, ( 1 ) 

where 7 ^ and uji {i = 1 , 2 ) are damping ratio and un¬ 
damped angular frequency of the i-th oscillator, and k 
is the coupling constant. Eigure[T] shows the time series 
of xi and X 2 of Eq. m in the logarithmic scale when 
uji = co ’2 = 1.0 and 71 = 0. Eirst, we consider uncou¬ 
pled case, k = 0. As we set 71 = 0 and 72 = 0 . 1 , the 
time series of xi exhibits a stationary oscillation without 
damping, while an exponential damping appears in the 
time series of X 2 , as shown in Eig. [IJa). Next, we con¬ 
sider a finite coupling strength of /c = 0 . 1 . In Eig. [T](b) 
with 7 i = 0 and 72 = 0 . 1 , both time series of xi and 
X 2 exhibit decays with envelope oscillations. Their de¬ 
cay rates, given by the slope of time series of xi and X 2 
in the logarithmic plot, are equal. As 72 increases from 
0 . 1 , the period of the envelope oscillation and the de¬ 
cay rate increase. At 72 ^ 0 . 2 , the envelope oscillation 
disappears and the decay rate reaches a maximum (see 
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FIG. 2 : (color online), (a) Real and (b) imaginary parts of 
two eigenvalues of which imaginary parts are positive as a 
function of 72 when uj 2 — 1.0 (black) and 0 J 2 — 1.005 (red) 
with 71 = 0.0 and cji = 1.0. (c) Real and (d) imaginary parts 
of the eigenvalues near EP as functions of ijJ 2 and 72 when 
uji = 1.0 and 71 = 0.0. The black circle, red dotted line, and 
blue dotted line represent the EP, real value crossing line, and 
imaginary value crossing line, respectively. 


The stability matrix M is then given by 


EIG. 1: (color online) Time series of xi (Black) and X 2 (Red) 
when cci = ct ;2 = 1.0 and 71 = 0.0. (a) No coupling case of 
/c = 0.0 with 72 = 0.1. Goupling cases of k = 0.1 with (b) 
72 = 0.1, (c) 72 = 0.2, and (d) 72 = 0.3. Insets show linearly 
scaled time series. 


Fig.[T]^c)). When 72 increases further, the decay rate de¬ 
creases again. For example, the time series of the case 
with 72 = 0.3 is shown in Fig. [1] (d). Although the am¬ 
plitude of two oscillators are different, as shown in the 
inset, their decay rates are equal. In our work, we con¬ 
centrate on the case that each uncoupled oscillators have 
zero or weak damping ratio so that their dampings are 
underdamped. 

In order to understand the variation of decay rate with 
72 and its maximum at 72 ^ 0.2, we analyze the eigen¬ 
values of a stability matrix around the origin. Eq. (HD 
can be rewritten as 

= yi, 

m = -712/1 - - kx2, 

X2 = 2 / 2 , 

m = -722/2 - w|a ;2 - kxi. (2) 

This set of equations is represented by a vector equation, 
z{t) = Mz(t), where z{t) = (xi{t),yi{t),X2{t),y2{t))'^■ 
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The eigenvalues Xi of M are complex numbers, because 
the matrix M is non-Hermitian. Since the time evolution 
of an eigenvector ei is given as ei{t) = ei exp(A/t), the real 
and imaginary parts of the eigenvalues correspond to the 
decay rates and the angular frequency of the correspond¬ 
ing time series, respectively. 

The complex eigenvalues with positive imaginary parts 
are shown as a function of 72 in Fig.[ 2 l^a) and (b). When 
72 < 0.2, real parts of two eigenvalues are very close but 
their imaginary parts are quite different, this means that 
the dynamics of eigenvectors would show almost same de¬ 
cay rate and different angular frequencies. In this range, 
the time series of xi and X2 would show a constant overall 
slope given by the close real parts, but they would have 
an oscillatory envelope whose frequency is determined by 
the difference of the imaginary parts of eigenvalues. This 
behavior has been shown in Fig. [TJb). As 72 approaches 
to 0.2, the real parts of two eigenvalues decrease and the 
imaginary parts become closer with each other, which 
corresponds to the time series with a faster decay and a 
longer period of envelope oscillation, respectively. 

As 72 goes further beyond 0.2, two real parts start to 
split but the difference of two imaginary parts become 
small. The splitting of two real parts indicates that the 
time series can be characterized by a combination of fast 
and slow decays. The fast decay might be seen only in 
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the short time behavior and the slow decay, correspond¬ 
ing to the larger real part, dominates the long time be¬ 
havior of the time series. Thus, although two imaginary 
parts are still different, there is no envelope oscillation 
due to the fast suppression of one eigen-component with 
the lower real part (see Fig. [U^d)). Note that the larger 
real part, governing long-time behavior, has a minimum 
value around at 72 ^ 0.2, which explains the maximum 
decay rate observed in Fig. [IJc). 

Note that two complex eigenvalues are very close at 
72 ^ 0.2 as shown by the black lines in Fig. [ 2 ] (a) and 
(b). We can expect that there should be a degenerate 
point, called exceptional point (EP) 0,0 , where two 
complex eigenvalues coalesce, in the system parameter 
space. By adjusting uj2 a bit as UJ2 = 1 . 005 , we find an 
EP at (co’2,72) ^ ( 1 . 005 , 0 . 2 ), which is shown by the red 
lines in Eig. [2] (a) and (b). It is well known that two 
eigenvectors also coalesce at the EP and mathematically 
the EP is the square-root branch point. The EP can be 
characterized by a peculiar eigenvalue surfaces in a pa¬ 
rameter plane. In Eig. [ 2 ] (c) and (d), the surfaces of the 
two eigenvalues are plotted in (002^12) plane. Topology of 
the surface explains the exchange of two eigenvalues for 
a parameter variation encircling the EP Q . It is empha¬ 
sized that the larger real part becomes a local minimum 
at the EP, indicating the local maximum decay rate in 
the parameter plane. 




FIG. 3: (color online). Maximal values of real parts of eigen¬ 
values with (a) R 2 = I.O, (b) R 2 = 0.0, and (c) R 2 = — I.O 
when Ri = 1.0. The colored and white region represent neg¬ 
ative and positive values, respectively. The blue dotted line 
represents the EP. (d) Real and (e) imaginary parts of two 
eigenvalues of which imaginary parts are positive as a func¬ 
tion of k when Au = 4.0 and Ri = 1.0. Black, red, and 
green curves represent the cases of R 2 = 1.0, 0.0, and —1.0, 
respectively. 


A. The amplitude death in coupled limit cycle 
oscillators 


III. EXCEPTIONAL POINT AND AMPLITUDE 
DEATH IN COUPLED LIMIT CYCLE 
OSCILLATORS 

In this section, we study the role of the EP when the 
amplitude death (AD) occurs in coupled limit cycle os¬ 
cillators. Let us start with the following system of two 
Stuart-Landau limit-cycle oscillators with diffusive cou¬ 
pling: 

zi = {Ri + iuJi — |2:i|^)zi T k{z2 — zi), 

Z2 = (i?2 + - 12:27^2 + k{zi - Z2), (4) 

where Zj are complex variables, ujj are the intrinsic an¬ 
gular frequencies of uncoupled j-th limit cycle oscilla¬ 
tors, and k is the coupling strength. Without coupling 
{k = 0), two limit cycle oscillators are attracted to the 
limit cycle with radii ^/^ for Rj > 0 and the origin for 
Rj < 0 . Stuart-Landau limit-cycle oscillator is renowned 
as a paradigmatic model for studying the AD in coupled 
nonlinear oscillators because it is a prototypical system 
exhibiting a Hopf bifurcation that can reveal universal 
features of many practical systems. Eor instance, a vari¬ 
ety of spatio-temporal periodic patterns can be created in 
two-dimensional lattice of delay-coupled Stuart-Landau 
oscillators 


It has been well known that the AD occurs in coupled 
limit cycle oscillators at proper k if the Aoo = 002—^1 is 
sufficiently large when Ri = R2 = I.O dMl. In order 
to obtain the AD region in the parameter space (Acj, /c), 
we calculate the Jacobian matrix J at the origin, which 
is given by 

/ Ri — k —uji k ^ \ 

j _ uji Ri — k 0 k I . . 

^ ~ k 0 R 2 -k -002 j ' ^ ^ 

\ 0 k UJ 2 R2 — k / 

The eigenvalues A of J are complex numbers because the 
Jacobian matrix J is a non-Hermitian matrix. That is, 
the real and imaginary parts are the decay (or growing) 
rates and the angular frequency of the orbit near the 
origin, respectively. 

The occurrence of AD is determined by the stability of 
the origin, which is related to the maximal value of the 
real parts of complex eigenvalues. If the maximal value is 
negative, the origin is stable fixed point and therefore the 
system exhibits the AD. The colored region in Eig. [ 3 l^a)- 
(c) where the maximal value is negative represent the 
AD regions when R2 = I.O, 0 . 0 , and —I.O, respectively, 
with Ri = I.O. As R2 decreases from I.O to —I.O, the 
AD region becomes larger. Eigure ^d) clearly shows 
the transition between positive and negative values of 
maximal real parts as a function of k when Auj is fixed. 
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FIG. 4: (color online), (a) Real and (b) imaginary parts of 
two eigenvalues near the EP at {AR,k) = (0.0,2.0) when 
Auj = 4.0 


B. The exceptional point in coupled limit cycle 
oscillators 

Similarly as the case of coupled damped oscillators, 
there also exists an EP in the coupled limit cycle os¬ 
cillators. The EP occurs at k = 2.0 when Auj = 4.0 
and Ri = R2 = 1.0, which is the double root position 
in Fig. [ 3 fd) and (e). Considering Ri = R2 = i?, four 
eigenvalues of Eq. m are given by 

-k + R± + 

-k + R± + + (6) 

where A = Aui^ — 4 fc 2 -(- 2 wi) +wi(\/ Aui^ — 4fc2 

wi) and □ = - 4k‘^-2aji)+ui{VAoj^ - 4P- 

cji), respectively. From the condition for EP, i.e., A = 
—the analytic condition for the existence of EP is 
given by 


Ri=R 2 , k = Auj/ 2 . ( 7 ) 

The eigenvectors also coalesce at this condition. Accord¬ 
ing to the Eq. ([ 7 ]), the EP occurs on the line in the param¬ 
eter space (Acj, k) when AR = R2 — Ri = 0.0 as shown 
in Fig. [ 3 l^a). If Auj is fixed, it is expected that a system 
shows the fastest attracting to the AD state on the con¬ 
dition of EP, k = Auj { 2 . Because the decaying rate to 
the AD state can be considered as a maximal value of 
Re(A) and the maximal value of Re(A) has its minimum 
at the condition of EP, k = Aco ’/2 [cf. Fig. [ 3 l^d)]. In 
addition, there is transition of transient behavior to the 
AD state on the EP, which is the transition between de¬ 
caying with envelope oscillation due to the effective beat 
note for k < Auj/ 2 and decaying without envelope oscil¬ 
lation for k > Auj/ 2 . It is noted that as R2 decreases 
from 1 . 0 , the AD region becomes larger, while the fastest 
attracting to AD state occurs on the EP when R2 = 1 . 0 . 

Figure m shows the complex eigenvalues near the EP in 
the parameter space (Ai?, /c), which is the singular point. 
We note that there are the same topological structures 
of eigenvalues near the EPs in coupled limit cycle oscil¬ 
lators with the parameter space {AR, Auj) when k = 2.0 
according to the Eq. O- 



EIG. 5: (color online). Time series of real parts of (black) 
and Z 2 (red) with (a) k = 0.5, (b) 1.1, (c) 2.0, (d) 2.4, and (e) 
3.0 when Auj = 4.0 and Ri = R 2 = 1.0. 


C. Numerical results 


In order to confirm the role of the EP expected in the 
previous subsection, we obtain the time series of zi and 
Z 2 as k increases. Figure [ 5 ] shows the time series of real 
parts of 2^1 and Z 2 with different k when Auj = 4.0 and 
Ri = R2 = 1-0. The coupling is turned on at t = 10.0, 
i.e., k = 0.0 when t < 10 . 0 . At k = 0 . 5 , neither AD 
nor 1 : 1 frequency locking occurs because of small cou¬ 
pling strength. Let us remind that the AD and 1:1 fre¬ 
quency locking occur, when the maximal value of Re (A) 
in Fig. [ 3 fd) is lower than zero and the pair of Im(A) in 
Fig. [ 3 fe) equal each other, respectively. At k = 1.1, the 
AD occurs with transient behavior of envelope oscillation 
but there is no frequency locking on the transient behav¬ 
ior. At /c = 2 . 0 , the AD occurs without envelope oscilla¬ 
tory transient behavior and the decay is fastest because 
this is the condition of the EP. The 1 : 1 frequency lock¬ 
ing on the transient behavior is also shown. At k = 2 . 4 , 
the AD occurs without envelope oscillation and there is 
frequency locking on the transient behavior. The decay 
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FIG. 6: (color online). —{I/Iad) as functions of Acj and k 
with (a) R 2 — 1.0, (b) 0.0, and (c) —1.0 when Ri — 1.0. 
The colored and white region represent the AD and non-AD 
regions, respectively. The blue dotted line represents the EP. 
(d) tAD as a function of k when R 2 = 1.0 (black), 0.0 (red), 
and —1.0 (green) when Acc = 4.0. 


is slower than that in the case of k = 2.0. At/c = 3 . 0 , 
the AD does not occur but there is frequency locking. 
In the AD region ( 1.0 < /c < 2 . 5 ), the EP is the tran¬ 
sition point between decaying with and without envelope 
oscillations. Also, in this region, the EP is the transi¬ 
tion point for frequency locking. The imaginary parts of 
eigenvalues relating to the frequencies change two differ¬ 
ent values into one value via the EP when R2 = 1 . 0 . If 
7^ ^2, two different frequencies are changed into two 
close frequencies not an identical frequency and therefore 
there is no exact frequency locking of transient. 

The important role of EP in AD is that the condition 
of EP guarantees the fastest attracting time to the origin, 
i.e., the AD state. We investigate the attracting time to 
the AD state, denoted by tAD- Here, tAD is calculated by 
followings; If, at time t, the radii of two oscillators firstly 
become smaller than cad^ i-e., the criterion for the AD 
state, and continue to be smaller than cad for 200 sec¬ 
onds, then tAD equals to t — ton where ton is the time 
when the coupling is turned on. We set cad = 0.001 and 
ton = 10 . 0 . Figure [6] (a)-(c) show —(l/t^i)), with vari¬ 
ous R2 when Ri = 1.0, on the parameter space (Acj, k). 
Figure [6] (d) shows — (I/^ad) as a function of k when 
Ri = 1.0 and Auj = 4.0 and the local minimum appears 
more clear when the parameters of system are closer to 
the EP. Contrary to the expectation from the maximal 
real parts of eigenvalues in Fig. [ 3 l there are many wrin¬ 
kled patterns when R2 = 1 . 0 . The wrinkled patterns 
gradually disappear as R2 decreases and then there is no 
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AR AR 

FIG. 7: (color online), (a) Maximum of real parts of eigen¬ 
values and (b) —{1/tAD) as functions of AR and k with 
Acc = 4.0. The colored and white region represent the AD 
and non-AD regions, respectively. 


patterns when R2 = — 1 . 0 . The different ton which means 
the different initial conditions makes the different wrin¬ 
kled patterns. The wrinkled patterns when k < Auj 2 
are caused by the oscillatory transient behavior. How¬ 
ever, the reason of the wrinkled patterns when k > Aco’/2 
is that the transition from fast decay to slow decay oc¬ 
curs when the amplitudes of the oscillators are smaller 
than our critical value cad and therefore the patterns 
disappear if cad is sufficiently small. 

Figure [ 7 | shows the maximal values of real parts of 
eigenvalues and —{ 1 /tAD) with the parameter space 
(Ai 7 , k) when Aoo = 4 . 0 . As shown in Fig. |H the EP 
exists at (Ai 7 , /c) = ( 0 . 0 , 2 . 0 ), where the maximal values 
of real parts of eigenvalues are local minimum as shown 
in Fig. [3 (a). tAD are also local minimum at the EP. In 
Fig. El^b), the wrinkled patterns exist at AR = 0.0 but 
they disappear as AR deviates from 0 . 0 . In principle, 
for the long time behavior, the oscillation behavior such 
as underdamped case exists only on the line, AR = 0 
and k < Aco’/2, because the real parts of two eigenvalues 
are same on the line in the parameter space (Ai 7 , k). 
Different real parts of two eigenvalues mean the system 
has two different decay rates and therefore only one fre¬ 
quency is dominant for a long time behavior. It is noted 
that the EP is not local minimal point on the parameter 
space {Auj, k) because the EP forms the lines as shown 
in Fig. [ 3 ] (a) and Fig. [6] (a). That is, the maximal values 
of real parts of eigenvalues decrease as the Auj increases 
on the EP line, k = Auj/ 2 . 

IV. SUMMARY 

We have studied the exceptional point in dynamical 
systems and investigated the role of the exceptional point 
in the transient behaviors of amplitude death in coupled 
limit cycle oscillators. The exceptional point is associ¬ 
ated with a critical point of frequency locking as well as 
the transition of the envelope oscillation, which also gives 
the fastest decay to the amplitude death in coupled limit 
cycle oscillators. In addition, for other examples (two 
Van der Pol oscillators interacting through mean-field dif- 
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fusive coupling, and coupled system of the Rossler and 
a linear oscillator), we have obtained the largest decay 
rates and transition behaviors at exceptional point (not 
shown here). As a result, the transient behaviors related 
to the exceptional point appear commonly for the cou¬ 
pled dissipative dynamical systems, independent of the 
specific properties of systems. We expect the exceptional 
point is important to the study on the various disciplines 
such as the noneqnilibrinm statistical mechanics and 
transient chaos [^, because the exceptional point is 
not related to the stationary states but the transient be¬ 
haviors. 
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